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We present results of the first experiment to probe for the electric dipole moment (EDM) of the electron using 
an £l-doublet state in a polar molecule. If the molecule is both massive and has a large molecular-fixed frame 
dipole moment, then the fl-doublet states have the potential to greatly increase the sensitivity of experiments 
searching for the EDM while also allowing for new methods of systematic error rejection. Here, we use the 
metastable a(l)''S+ state of lead monoxide (PbO) to probe for the electron EDM. Our best fit for the electron 
EDM of = (-4.4 ± 9.5stat ± l-8syst) x 10"^' e-cm allows us to place an upper limit on the magnitude of the 
EDM of \di,\ < 1.7 X 10"^* e-cm (90% confidence). While this is less stringent than limits from other, previous 
experiments, our work emphasizes the systematic error rejection properties associated with the £l-doublet level 
structure. The results should inform the work of other, ongoing experiments that use molecules with analogous 
level structure. 

PACS numbers: 14.60.Cd, 33.40.+f, 11.30.Er, 33.57+c 



I. INTRODUCTION 

The detection of a permanent electric dipole moment 
(EDM) of a fundamental particle would provide evidence of 
violation of charge-parity (CP) symmetry [1]. If an EDM ex- 
ists, it violates both parity (P) and time (T) reversal symme- 
tries. Through the CPT theorem, which states that the com- 
bined operations of parity, time, and charge-conjugation (C) 
must be conserved in any Lorentz-invariant theory [2], a vio- 
lation of T symmetry, such as that by an EDM, is equivalent 
to violation of CP symmetry. CP violation was first observed 
in the decay of the neutral kaon [3], and such violation can 
be explained through the standard model (SM). However, the 
SM does not contain enough CP violation to explain the cur- 
rent matter-antimatter asymmetry in the universe [4]. Theo- 
ries that go beyond the standard model generally provide for 
more CP violation, and therefore, larger EDMs [5-9]. 

No EDM of a fundamental particle has yet to be detected; 
however, a recent experiment using YbF molecules has set 
the current experimental limit on the electron EDM (eEDM) 
of \dA < 1.05 X 10-2^ e-cm (90% confidence) [10]. This ex- 
periment, like most of its predecessors [11, 12], attempted to 
detect the eEDM by detecting a change in Larmor precession 
in the presence of an electric field. To understand how these 
experiments work, consider a free electron in the presence of 
both electric and magnetic fields. If the electron possesses an 
EDM, the magnitude of the resulting Larmor precession fre- 
quency vector u) is given by 

CO ^ j\2g,fiBR + 2de£\ , (1) 



where £ is the applied electric field, dg is the eEDM, /ig is 
the Bohr magneton, ge is the g-factor for the electron, B is the 
magnetic field, and H is Planck's constant divided by 2:7r'. If 
the two fields are nominally applied along the same axis and if 
the EDM is parallel to the magnetic dipole moment^, the mag- 
nitude of the precession frequency will be larger (smaller) if 
the two fields are parallel (anti-parallel). To isolate the preces- 
sion due to deS, one typically measures the difference in the 
precession frequency when the relative direction of the S and 
B fields are reversed. 

Because the electron is charged, it is not feasible to mea- 
sure the precession frequency of free electrons, and therefore 
atoms and molecules with unpaired electron spins are com- 
monly used [13, 14]. If an electric field is applied to such 
an atom or molecule, Eq. 1 can still be applied, except 
is replaced by the ^-factor of the bound state with which the 
measurement is performed, and S is replaced by an effective 
electric field Seff with which the electron spin interacts. In 
the non-relativistic limit, Seff is equal to zero (a result known 
as the Schiff theorem), but due to relativistic effects can be 
non-zero or even larger than the applied electric field. Atoms 
or molecules with a heavy nucleus can offer effective electric 
fields as high as 100 GV/cm. Because de < 10"^^ e-cm, the 
shift in Larmor precession frequency will still be small, of or- 
der 2;r X 0.1 Hz. The bias magnetic field B must still be large 
enough to observe precession; in practice, B is ordinarily of 
the order of 10 to 10^ mG, corresponding to overall precession 
frequencies in the range of w ~ 27r x (10"*) to 27r x (10^) Hz. 
Therefore, the eEDM may produce a tiny fractional shift of 
between 10""^ and 10"** in the spin precession. 

In any experiment where the expected signal represents 
such a small fractional shift, the experimentalist must take 
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We note that the precession frequency is a vector in the most general case. 
While the experiment discussed in this work only measures the magnitude 
of the precession frequency vector, there are experiments that can measure 
not just the magnitude, but the direction of the precession. The implications 
of such a measurement are discussed briefly in Appendix A. 
- I.e., both are anti-parallel to the spin. Such a condition would imply < 0, 
just as the gyromagnetic ratio = gefiB obeys < 0. 
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care to isolate systematic effects that mimic the sought-for 
shift. For eEDM experiments, a magnetic field that is gen- 
erated by the electric field, such as one created by leakage 
current that flows from one electrode to another, represents 
one of these possible systematic eff'ects. 

In experiments where the Larmor precession frequency is 
given by Eq. 1, reversal of the magnetic and electric fields is 
the primary method used to measure and reject systematic ef- 
fects. In this method, magnetic and electric fields are applied 
along a chosen laboratory axis and their directions relative to 
that axis are reversed. Four Larmor precession frequencies 
corresponding to the four possible applied field configurations 
are then measured. Combinations of these measured frequen- 
cies then yield information regarding the eEDM and the ap- 
plied fields. For example, if ojij represents the measured pre- 
cession frequency when B has sign / and & has sign j, the sum 
of aU four obeys dLi+ + + + + w_ _ - ^gfisBav, where 
Bav is the average magnetic field magnitude. The eEDM can 
be found by forming the combination (y+ + -co+- - H-a;_,_, 
which assuming all the field magnitudes are identical, yields 
%de&eff- Given that there are four different measured fre- 
quencies, there are two other field parameters, i.e., parame- 
ters that specify the applied fields. These correspond to a non- 
reversing component of the magnetic field B„r and a change in 
the magnetic field magnitude Bcon- that is correlated with the 
absolute direction of the electric field. However, a magnetic 
field component whose sign depends on both the sign of S and 
B, e.g., a magnetic field that is generated by leakage current 
that flows between the two electrodes, will be completely in- 
distinguishable from eEDM in this method. In addition, no 
information regarding the electric field or the alignment of the 
two fields can be obtained solely from the measured preces- 
sion frequencies. 

Recently, considerable attention has turned to polar 
molecules with Q-doublet substructure [15-19], because such 
systems have the possibility to reject systematics better 
than the simple atomic or molecular experiments described 
above [20]. Such O-doublet substructure occurs generically in 
molecular states with internal electronic angular momentum 
Je > 1 • In addition to having states with Q-doublet substruc- 
ture, the molecular species must generally have one heavy nu- 
cleus, two valence electrons in a triplet state, and at least one 
valence electron in a cr or tti /2 orbital in order that Se/f for the 
Q-doublet states is also large. As shown in Fig. 1, it is possi- 
ble to prepare such a molecule in quantum states that have op- 
posite signs of the molecular polarization and hence opposite 
signs of the effective electric field. This additional degree of 
freedom can possibly yield more information regarding exper- 
imental conditions such as field alignments, field magnitudes 
and leakage currents'*. 

A detailed description of the level structure of these states 
is given in Refs. [24, 25], and a short review will be presented 
here. With no applied magnetic and electric fields, the total 



^ Other potential systems are being studied for eEDM searches as well, in- 
cluding solid state systems [21, 22] and Cs and Rb atoms in an optical 
lattice [23]. 



M = -l M = M = +l 




FIG. 1. (Color online) Schematic of the energy level structure of 
a y = 1 f2-doublet in a polar molecule. (Top panel) With no ap- 
plied electric or magnetic fields, there are two sets of three degen- 
erate states, coiTesponding to the parity eigenstates P and the three 
projections M of the angular momentum J along the quantization 
axis, respectively. Which of the parity eigenstates has higher energy 
depends on the exact structure of the molecule in question; in the 
a(l)'E^ state of PbO used in this work, the P = +1 state lies lower 
in energy. (Bottom panel) An applied electric field £ (red, long dash 
vector) mixes the opposite-parity states (shown by the dotted, gray 
lines) and induces a Stark shift/splitting between them. The Stark- 
shifted levels (shown by the red, long dashed lines) have a net aver- 
age molecular dipole moment n = ji„h, where h is the molecular axis, 
oriented either along or against the electric field. Hence the field £ 
naturally defines the quantization axis. An applied magnetic field B 
(blue, dash-dot vector) in the same direction as £ induces a Zeeman 
shift between the M = ±1 states (shown by the blue, dash-dot lines). 
The eEDM interacts with the effective electric field of the molecule, 
which has opposite sign for the two A' states, and modifies the Zee- 
man shift (black lines). Last, the electric field causes mixing between 
the J = I and J = 2 states that effectively changes the g factors of 
the upper and lower A'-states by an amount \k\S\. 

angular momentum 7=1 state has two degenerate manifolds 
of three states. Each of the three states corresponds to a dif- 
ferent projection M of the total angular momentum along the 
laboratory z axis, as shown in the top panel of Fig. L The 
two sets of degenerate states correspond to eigenstates"* of the 



* These parity states are best expressed in the basis \J, M, \Cl\ = 1, P = ±1), 
where il is the projection of the electronic angular momentum Jg along 
the molecular axis n, i.e., £1 = Jg ■ h. The paiity states can be expressed 
in the signed n basis through \J. M, |n| = 1, P = ±1> = -L(|7, M, £1 = 1> - 

i-iyP\J,M,a. = -1)). 
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parity operator with quantum number P - ±1. Because of 
a Coriolis coupHng to the total electronic angular momentum 
Jc, these two parity states have a small difference of energy 
An. Which of the two parity eigenstates has higher energy de- 
pends on the specific molecule and state. For the purposes of 
this discussion, the state with higher energy will be denoted 
by H, and the state with lower energy will be denoted by L. 
Note that while in the molecule-fixed frame there is an elec- 
tric dipole moment ju^, the expectation value in the laboratory 
frame for either of these states is identically zero. In addition, 
the H and L states have slightly different magnetic ^-factors, 
which will be denoted as gn and gi, respectively. 

An applied electric field causes the molecules to polar- 
ize along or against the electric field. The tensor Stark shift 
that accompanies this polarization is shown schematically in 
Fig. 1. In the presence of £, the projection of the molecular- 
fixed dipole moment fig along the electric field, given by 

= sign{{ fia) ■ £), becomes a good quantum number Thus, 
the direction of the electric field defines the quantization axis 
for the molecule^. When a magnetic field is applied along 
the quantization axis, it interacts with the unpaired spins and 
causes a splitting between the M = ±1 sublevels of IgN^jUsB, 
where gN+{N-) is the g-factor for the A^=1(A^ = -1) state. As 
described in Sec. VII of Ref. [24], the g factors diverge when 
strongly polarized by an applied electric field, i.e.. 



M - +1 states is given by 



20Br 



(gH + gl) 



(2) 



where B,. is the rotational constant. 

The eEDM itself interacts with the effective electric field of 
the molecule, not the external field. Because the states have 
different signs of A^, we write the energy shift between the 
M - +\ states due to the eEDM as £ = -2de ■ where 
E>eff - —£>effn is the effective electric field with which the 
eEDM interacts; this quantity has opposite sign for the two A^ 
states. Once the applied electric field is sufficiently large to 
polarize the molecule, fi^// reaches its maximum value and 
becomes independent of the applied electric field^. 

Because a polarized molecule in the J - \ state has a large 
energy difference between the M - Q and M - +\ states, 
magnetic fields perpendicular to the quantization axis, defined 
by the electric field £, will have a minimal impact on the pre- 
cession frequency. Such a transverse magnetic field Bj will 
couple states with M and M' - M + 1, and therefore there 
is no shift to first order in Bj. To second order, assuming a 
fully polarized molecule and neglecting the effect of the J - 2 
states, the shift in energy 6E due to Bj for the four A^ = ±1, 



useful to 
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^ For this configuration, it is 
|7, M, N) = \j. M,a = NsigniM)) 

(-l)"|/,M,|n| = 1,P = 
* In the case of a paramagnetic atom, G^ff is hnearly proportional to the ap- 
plied electric field for all fields achievable in the laboratory. This dilference 
between atoms and molecules is notable and could possibly be used as an 
additional check for systematic eft'ects. 
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where B; is the magnetic field along the electric field, Bt is the 
transverse magnetic field perpendicular to it and g = ^{gH + 
gl). Therefore, while a transverse field may lead to a shift 
in the overall energy of the M = ±1 states, such a shift is 
mostly common mode. The actual change in the precession 
frequency for the state with quantum number A^ due to the 
transverse magnetic field dcj^ is given by 



J. dEN.M=+i 
5a>j^ = 



^-2{gH + gl)- 



(4) 

which is the same for both A^ = +1. 

Collecting all the terms above, we can write the magnitude 
of the precession frequency of the A^ = ±1 states to lowest 
order in the various fields as 
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+ 2Nde6eff 



(5) 



where k - 3jj.a/lOBr. Unlike experiments based upon Eq. 1, 
it is now possible to measure eight different Larmor preces- 
sion frequencies corresponding to the four different applied 
field configurations and two different N states. These eight 
measured frequencies are denoted by ojij^k, where ; = sgn(A^), 
j = sgn(£) and k = sgn(B), and the eEDM can be found 
by taking the combination of the measured frequencies that is 
odd in A^, £, and B, i.e., ij^'^ij.k - l6dgSi,ff- 

Moreover, an experiment based on measuring the eight dif- 
ferent Larmor frequencies given by Eq. 5 can yield more infor- 
mation regarding the applied fields. First, with applied elec- 
tric fields S < 100 V/cm, kS is typically of the order of 10"^. 
Therefore, the two A^ states respond almost equally to the ap- 
plied magnetic field and act as an internal comagnetometer, 
allowing cancellation of the effects of a fluctuating or system- 
atically changing magnetic field. For example, this can be 
used to distinguish frequency shifts due to a field generated by 
current flowing between the electrodes from shifts due to the 
eEDM. Second, while the electric field dependence of the g 
factors prevents perfect cancellation of magnetic field effects, 
it does allow extraction of information regarding the electric 
field. This information can be used to determine whether the 
electric field magnitude is changing with time or upon reversal 
of its direction. 

An experiment based upon Eq. 5 also can reject systematic 
errors due to misalignment of the fields. With typical align- 
ment errors between the electric and magnetic fields, the ratio 
jU^B x£/(//fl£) can be of the order of 10"^, leading to rejection 



of transverse magnetic fields of the order of lO"'". Moreover, 
motional magnetic fields generated by the relativistic v x £ 
effect will have a negligibly small impact on the eEDM mea- 
surement. 

In this work, we examine the systematic rejection prop- 
erties of an eEDM search based on this level structure, us- 
ing an experiment based on the a^I,^ state of lead monox- 
ide (PbO). PbO has proven interesting for an eEDM search 
for a multitude of reasons. First, it can be produced with 
high-densities in a vapor cell [24, 26-29]. Second, it has 
both large electric dipole moment {fia/fi ~ In y. 1.64 MHz 
V^'cm"' [27]) and a large effective electric field (£f// ~ 
25 GV/cm [28, 30]). While the state is sensitive to magnetic 
fields, with g ^ 0.86 [29], the difference in the g factors with 
magnetic fields obeys k ^ 1 x 10"^ V^'cm"', where we have 
used Br in = 27r X (7.054 GHz) [31]. This value of k ensures 
that the molecule will behave as a good comagnetometer for 
our values of applied electric fields. 

Although the present work was not able to establish a new 
limit on the eEDM, many future experiments depend on the 
capabilities for the systematic rejection of the level structure 
shown in Fig. 1, and we believe the findings of this work will 
prove useful for those experiments. We first discuss the appa- 
ratus' used in this study in Sec. II. The suppression of effects 
due to transverse fields is shown experimentally in Sec. III. A 
comprehensive list of frequencies or phases measured in such 
an experiment, along with a formal definition of the measur- 
able systematics and physical quantities, will be developed in 
Sec. IV. In addition to the effects detailed in Sec. IV, we con- 
sider the effect of field gradients on the experiment in Sec. V; 
these prove to be the dominant source of systematic errors in 
this work. The final result of the experiment is detailed in 
Sec. VI. Note that for the remainder of this article, we shall 
take h.= \ and treat energy and angular frequency as identical 
quantities. 



II. APPARATUS 

For our experiments, a gas of PbO molecules is created in 
a high-temperature vapor cell, as shown in Fig. 2. The cell 
is shaped roughly like a cube with approximately 3 in. sides 
and is made from fused alumina, which forms the walls and 
structural supports. Four, 2 in. -diameter, yttrium-aluminum- 
garnet (YAG) windows provide optical access from four sides 
of the cell. Re-entrant electrodes, with an approximate diame- 
ter of 2.25 in., protrude into the top and bottom with a spacing 
of 1.5 in. These electrodes are made from gold foil which 
has been adhered onto a flat, sapphire substrate using oxide 
bonding [32]. Surrounding the main electrodes are guard ring 
electrodes of inner diameter 2.5 in. and outer diameter 2.75 in. 
When a larger voltage is applied to these guard rings than the 
main electrodes, the electric field becomes more uniform in 
the main volume of the cell. 
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' For a comprehensive description of tlie appartus and basic experimental 
procedure used, see also Ref. [25]. 



FIG. 2. (Color online) A schematic of the experiment. The cell is 
placed at the center of vacuum chamber (VC). Light from a laser en- 
ters through a quartz light-pipe to illuminate the molecules, and the 
resulting quantum beat fluorescence is observed using photomulti- 
plier tubes (PMTs). Microwaves can also be applied to the molecules 
through a microwave horn. The vacuum chamber and cell are con- 
tained within three layers of magnetic shielding. The magnetic field 
coils and the oven that heats the vapor cell are not shown. The green, 
dotted arrows indicate the direction of rotation of the resulting quan- 
tum beat fluorescence. 



An oven heats the vapor cell to ~ 700°C, where the par- 
tial pressure of PbO is ~ 10 Torr and the total pressure 
of PbO vapor (dominated either by Pb404 or PbO) is ~ 
10"^ Torr [33, 34]. The heating elements are formed from 
laser-cut tantalum foil, designed to minimize the overall self- 
inductance, and are held in place by a quartz structure. Be- 
cause tantalum will oxidize when heated in atmosphere, the 
oven and vapor cell are placed within an aluminum vacuum 
chamber with 18 in. height and 14 in. diameter, where pres- 
sures of approximately 10"^ Torr are achieved. Two-inch di- 
ameter quartz light pipes protrude through the walls of the 
vacuum chamber and provide the necessary optical access for 
excitation and detection of the molecules. 

The vacuum chamber is surrounded by multiple magnetic 
field coils. A 300 turn, 10.5 in. radius Helmholtz coil gen- 
erates a magnetic field in the z direction. Cosine-type coils 
apply uniform fields in the x direction and y direction. A set 
of gradient coils can generate all possible linear magnetic field 
gradients. The vacuum chamber and magnetic field coils are 
placed within three layers of /i-metal shielding, which provide 
a shielding factor of the order of 10^. 

A short laser pulse of 548 nm light prepares the molecules 
into the a^Z^ state of PbO. The pulse is generated by a nar- 
rowband, continuous-wave-seeded, pulse-pumped dye ampli- 
fier [35, 36]. A diode laser, which is amplified using a semi- 
conductor tapered amplifier and frequency-doubled using a 
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periodically-poled lithium niobate (PPLN) waveguide, gen- 
erates the seed laser light. An injection-seeded, frequency- 
doubled Nd: YAG laser with an output pulse of approximately 
3 ns in duration and a repetition rate of 100 Hz acts as the 
pump for Fluoresin 548 dye. The output energy of the full dye 
amplifier system has an energy of approximately 7-10 mJ per 
pulse. Because the dye amplifier's linewidth (< 200 MHz) is 
smaller than the Doppler width of the molecular transition at 
T K 1000 K (Td ~ 800 MHz), a mirror retroreflects the beam 
to excite molecules in a broader range of velocity classes. 

With i-polarized light, the laser drives the transition from 
the absolute ground state X'2(v = 0, 7 = 0) into the mani- 
fold of ^E+(v - 5, J - 1) sublevels and coherently populates 
the = +1, M - +1 levels shown in Fig. 1. Neither the 
M = + 1 Zeeman splitting nor the = ± 1 Stark splitting is 
resolved within the Doppler width of the transition. Because 
of the non-negligible electric field inhomogeneity (~ 1%), any 
coherence between the = +1 states is quickly lost^. The 
resulting fluorescence signal indicates that the a state has an 
effective lifetime of approximately t„ x 50 fis (see below for 
details), with decay presumed to be due in nearly equal parts 
to spontaneous emission and quenching on cell walls. As the 
molecules decay to the ground state, quantum beats are cre- 
ated by interference of the decay paths to M = ground 
state sublevels [24]. The modulated fluorescence signal is de- 
tected using two photomultiplier tubes (PMTs) that are ori- 
ented along the i-axis, which is perpendicular to the laser's 
direction of propagation. Two filters mounted in front of each 
PMT, a KG-4 infrared-blocking colored glass and a custom 
554 +104 nm optical interference filter, serve to block the pri- 
mary spectrum of the substantial blackbody radiation present 
at r » 1000 K. The quantum beat fluorescence along the x- 
direction is polarized along the y axis; therefore, wire grid 
polarizers are also used to filter out other background light. 

In order to resolve the signal from the two different A' 
states, one of two techniques is used. As described in Sec. 
X of Ref. [24], microwaves that are resonant with either 
the 1/ = 1,M = +1,A^ = -1) -> |7 = 2,M = ±2,A^ = -1) or 
|7 = 1,M = +1,A^ = +1) ^ |7 = 2,M = ±2,A^ = +1) are ap- 
plied to the molecules for a duration tm » I/Qr, where Qr 
is the Rabi frequency of the microwave drive. This causes de- 
coherence of the beat signal from the particular A^ = + 1 state 
with which the microwaves are resonant. The frequencies of 
these transitions are approximately 28 GHz. The microwaves 
are generated using a custom-built microwave source, and are 
applied to the molecules via a microwave horn. 

Alternatively, one can also resolve the signals of the two A^ 
states by applying large magnetic and electric fields. If both 
fields are sufficiently large, the difference in the precession 
frequency between the two A^ states can become larger than 
l/T^, where is the lifetime of the quantum beats. In this 
case, both A^ states contribute to the quantum beat signal, and 
the two frequency components can be resolved, as shown in 
Fig. 3. In this case, there is no loss of signal due to microwave 
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FIG. 3. (Color online) An example of the resulting two beat signal 
from a single laser shot. At sufficiently large applied magnetic and 
electric fields, the difference in the g factors coupled with the size of 
the average beat frequency allows for resolution of both Q-doublet 
states at the same time, (a) The two precession frequencies generate 
a beat note in the raw signal (shown as red points), as evidenced by 
the node near f = 38 /js. This data fits well to Eq. 7 (shown as a blue 
curve), (b) At the beginning of the decay, oscillations at the average 
of the two frequencies are seen, (c) The beat signal still fits to Eq. 7 
after the node, where the phase of the signal has shifted by n. (d) 
A discrete Fourier transform with a rectangular window function of 
both the data (shown as red points) and the fit (shown as a blue curve) 
shows two resolved frequencies in the spectrum of the quantum beat 
signal. 



preparation, and therefore the sensitivity of the experiment to 
the eEDM is generally better in this two-beat case than in the 
microwave-erasure technique described above. 
The resulting quantum beats are fit to the function 



S(t) - b(t) + e a\ cos(wif -i- ^i) 



(6) 



or 



S(t) - b{t) + e [fli cos(wif + 4>\) + fl2Cos(w2f + ^2)] , (7) 

for microwave-erasure and two-beat techniques, respectively. 
In these fits, a,, 0,, w, and T are free parameters, f is the 
time since the laser pulse, and b{t) is the background signal. 
This background signal is determined directly from the data 
by using a low-pass, zero-phase-shift filter (a digital filter that 
moves both forward and backward through time [37]). In gen- 
eral, b{t) can be fit to 



b(t) = JlAe'"'" + Jlae'"'" + b 



•bh ■ 



(8) 



' We estimate that the N = ±1 coherence is lost in less than 250 ns. 



where t, and bhb are tunable constants. The constant term 
bht, allows for fitting the blackbody background. The first. 
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FIG. 4. (Color online) Evidence of electron emission. The current 
flowing into/out of the top electrode (TE) and its guard ring electrode 
(TGR) is shown as a function of the voltage applied to the TE, while 
the TGR is held at ground. When the applied voltage is greater than 
zero, the current flow shows an approximately Ohmic-behavior, with 
current from the TE being detected flowing into the TGR such that 
the sum is zero. When the voltage is negative, however, the sum is 
not zero and an excess of negative current is observed flowing from 
the top electrode. Directly above the cell is a heating element that is 
not in physical contact with the vapor cell. By biasing this heating 
element to -100 V (solid lines) or V (dashed lines), the threshold 
voltage at which excess current is observed is changed. This process 
is analogous to the change in current flow provided by an electron- 
emitting filament and a biased grid electrode in a vacuum-tube triode. 



fast exponential term generally obeys ta ~ (l/F) x 10" and 
is believed to be fluorescence from molecules that were ex- 
cited to the nearby A^n electronic state. The second, slower 
exponential decay obeys t„ < 1 /F and is believed to be the 
fluorescence lifetime of molecules excited to the a state that 
are not spin-polarized in the J=l state and hence do not con- 
tribute to the quantum beat signal. This fluorescence lifetime 
is the aforementioned x 50 fis. 

Using a vapor cell at 700°C poses some unique challenges. 
First, in order to maintain the temperature, current must flow 
through the heating elements. Such current will generate mag- 
netic fields that could interfere with the eEDM measurement. 
For this reason, the cuiTent through the heaters is brought to 
zero before the laser pulse is fired and subsequently restored 
approximately 400 //s after the laser pulse. In order to avoid 
inducing eddy currents in the aluminum vacuum chamber, the 
heaters are supplied with a 10 kHz AC signal, which during 
1 ms-long turn-on and turn-off periods is modulated by an en- 
velope of the form 1 - sin^(^f), where ^ = 2nx(2 kHz) [38]. 
Stereo audio amplifiers are driven with an arbitrary function 
generator to provide this AC signal with the ~ 1.1 kW of 
power necessary to heat the vapor cell. 

Secondly, operating at 700°C greatly decreases the electri- 
cal resistivity of most insulators, including the fused alumina 
and YAG used in our vapor cell. While this lower resistivity 
generally creates large leakage currents, an even more insid- 
ious problem plagues the vapor cell. At such large tempera- 
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FIG. 5. (Color online) Measurement of the alignment of the elec- 
tric and magnetic fields. In this polar projection plot, the radial and 
azimuthal components correspond to the polar angle and azimuthal 
angle, respectively, in a spherical coordinate system where +z is de- 
fined as the nominal direction of the magnetic field. With a fixed 
electric field applied to the molecules, the applied magnetic field 
was tilted from the J axis without changing its magnitude, and the 
resulting average beat frequency a>ay was measured at the colored 
points shown. The circular contours represent the best fit 2-D poly- 
nomial to the data. The place where the beat frequency maximizes, 
at 6) = (0.64 ± 0.02)° and ip = (128.2 ± 0.1)°, represents the place 
of maximal alignment between the magnetic and electric fields. The 
fractional uncertainty in the beat frequency measurement is approxi- 
mately 10"', using 512 laser shots per point. 



tures, electron emission is observed, as shown in Fig. 4. This 
electron emission and associated voltage drops across the path 
between the electrode leads and the emission surface have the 
possibility to seriously distort the electric field, e.g., changing 
its magnitude and direction. Making eEDM measurements in 
the presence of this problem provides one of the most rigor- 
ous tests of systematic quantification and rejection, and thus 
is integral to the primary goal of this work. 



III. REJECTION OF TRANSVERSE MAGNETIC FIELDS 

According to Eq. 5, the precession frequency is determined 
primarily by the projection of B onto the electric field £, 
which determines the quantization axis of the molecule. Any 
transverse field only affects the precession frequency at higher 
order For a typical experiment, the condition jUb(B x S) 
fia£' is satisfied, and the relevant transverse-field term in Eq. 5 
can be neglected. The precession frequency, averaged be- 
tween both = + 1 states, can then be written as 

COav = ^(^^+ = 2^//b(B • £) , (9) 
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where dLt± is the precession frequency for the = ±1 state^. 

If the alignment between the magnetic and electric fields is 
changed, Eq. 9 implies that the average precession will change 
as well, even if the magnitude of B remains constant. Such an 
effect is shown in Fig. 5. Here, an electric field is applied 
to polarize the molecules. The magnetic field is then tilted 
from its nominal direction without changing its magnitude. 
(In this coordinate system, the vertical or z axis is defined by 
the magnetic field applied by the Helmholtz coil, and x is the 
direction defined by the magnetic field generated by the cosine 
coil most aligned with the axis of the detectors.) If the angle 
between the two fields is denoted by y, then Eq. 9 becomes 

a>av - Igi^BB cosy - 2giJ.BB(cos 6 cos 6e 

+ sin0sin(9£ cos(0 - , (10) 

where 0e and (pE are the polar and azimuthal angles of the 
electric field, respectively, and 6 and cf) are the polar and az- 
imuthal angles of the magnetic field, respectively. If both po- 
lar angles are small, this can be expanded to second order and 
compared to a second-order polynomial fit to the data to ex- 
tract the relative angles. The best fit polynomial is shown as 
the contours in Fig. 5, and is accurate to within the errors of 
our calibration of the magnetic field. 

The fit in Fig. 5 not only justifies the approximation that 
misaligned components of the magnetic field are small com- 
pared to the primary component, but also demonstrates that 
the relative angles of the electric and magnetic fields in the 
experiment can be measured with much better than 1 ° preci- 
sion. 



IV. FIELD REVERSALS AND FREQUENCY 
COMBINATIONS 

As discussed in Sec. I, there are eight diff'erent precession 
frequencies W/j j that can be measured in a molecule with D.- 
doublet substructure. The indices / = sgn(A^), j - sgn(£) and 
k - sgn(B) specify the A' state measured and the field config- 
uration relative to a chosen laboratory axis. As in the case of 
an experiment based on Eq. 1, these eight different frequencies 
can be summed together to form a frequency combination. Let 
us define a combination of measured frequencies cjij^k as 

A/j-.A-w = (^+.+.+ + i(^-,+.+ + joJ+-,+ + ku+^+- (11) 
+ij<jj--^+ + ikto-,+ - + jka>+-- + ijka>--- . 

Note that because the frequency shift due to the eEDM is odd 
in all three reversals [see Eq. 5], A_,_,_w will correspond to 
the eEDM signal channel. Also, while this discussion focuses 
on frequency measurements, the following discussion applies 
equally well to experiments that measure the phase preces- 
sion of molecules in a beam. The differences between the two 
types of experiments are detailed in Appendix A. 



A. N-even combinations 

To determine the physical meaning of each combination, 
let us first consider the combinations that are even under re- 
versal of A^: A+_+_+w, A+___+aj, A+_+__aj, A+____aj. Because the 
g-factors of the = +1 states are affected by the electric field 
in the equal but opposite way, these combinations are nomi- 
nally independent of the electric field applied. For this reason, 
they give information only regarding the magnetic field. 

There are four total electric and magnetic field configura- 
tions and four combinations that are even in A^. Let us ex- 
amine each combination in turn. First, the term A+ + +w is 
proportional to the precession frequency averaged over both 
Q-doublets and the four states of £ and B. Therefore, we 
define A+ +_+w = 16gfiBav, where Bai, > is the average 
magnitude of the magnetic field B projected onto £. A non- 
zero A+.+._(y, corresponding to a shift of the precession fre- 
quency upon reversal of the magnetic field, can be due to a 
non-reversing magnetic field B„, , i.e., A+ = l6gfiBBnr- It 
is also possible that the electric field can produce a magnetic 
field whose direction is correlated with the direction of £ that 
will add vectorially to the applied magnetic field. Depend- 
ing on the direction of the applied magnetic field, this addi- 
tional field may lead to an increase or a decrease in the overall 
magnitude of the magnetic field. The average precession fre- 
quency will then change when both the electric and magnetic 
fields are reversed. The most obvious source of this additional 
£-induced magnetic field is due to leakage currents; there- 
fore, the associated frequency combination shall be denoted 
as A+ = \6giiBBieak- There is one additional combina- 
tion, namely A+ _ +w. A non-zero A+__ +w could arise if there 
is a change in the magnetic field magnitude that is directly 
correlated with the sign of the electric field, but not with the 
sign of magnetic field'°. The associated frequency combina- 
tion shall be denoted by A+ _ + - IGgfisBcorr- 

Given the above field parameterization, we can now ask 
what is the measured magnetic field magnitude for a given 
field configuration? For a particular field configuration, the 
measured, average frequencies between the A^ = ±1 states 
{(^av.j.k - \ Tji '^i.j,k) can be used to define the measured 
magnetic field magnitude B„,eas,i.k through the relationship 
<^av,iM - "^-gfJ^BBmeas.j.k- The relationships between the four 
magnetic field parameters {Bav, Bcrr, B„,, and Bieak) and 
Bmeas,j,k are shown in Table I. 

Here we have parameterized the frequency combinations 
in terms of magnetic field components often discussed in the 
context of additional fields in a traditional eEDM experiment. 
However, there is another way to understand what might cause 
non-zero values of A+ _ A+ + _<y or A+ _ As discussed 
in Sec. Ill, a change in precession frequency can occur if the 
axis of the magnetic field forms an angle with the electric 
field and this angle shifts upon reversal of either field. To 
understand how the relative angle might generate a non-zero 



' For the present discussion, we do not reverse the magnetic or electric fields, 
so we omit the j and k indices present in Section I. 



This could arise, for example, due to unwanted electronic coupling be- 
tween the magnetic and electric field power supplies, e.g., coupling due to 
a ground loop. 
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FIG. 6. (Color online) The angle interpretation for S,„, Bcorr, and 
Bleak- The electric field direction is defined by two vectors £^ and 
£l that correspond to it being in the positive and negative directions, 
respectively. The magnetic field is defined by two vectors (arbi- 
trarily defined to be the z direction in this plot) and that correspond 
to it being in the positive and negative directions, respectively. If no 
direction is equal, there are four angles formed between the electric 
field axes and magnetic field axes. These angles are specified by Oj^, 
where / and k correspond to the nominal sign of £ and B, respec- 
tively, in the laboratory frame. 
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TABLE I. The parametrization of the measured magnetic field mag- 
nitude B„,eas and measured electric field magnitude Gmeas, for the four 
combinations of the applied fields. The fields are nominally applied 
along a chosen laboratory axis, and j = sgn(£) and k = sgn(5) spec- 
ify the direction relative to that axis. For both the magnetic and elec- 
tric fields, the average fields (S„,, and £„,,) are assumed to be both 
larger than zero and much larger than the magnitude of any other 
component. Note that for the case of the electric field, there are three 
parameters and four field configurations. The missing parameter has 
the same symmetry properties as the eEDM, and is therefore indis- 
tinguishable from it. 



A+__ +w, A+ or A+____w, consider the specific case shown 
in Fig. 6. Let the magnetic field in the positive (negative) 
state be defined as the vector Bkd and the direction of the 
electric field in its positive (negative) state be defined by the 
vector When both the magnetic and electric fields are 

in the positive state, the measured precession frequency aver- 
aged over the = ± 1 states will be given by 

= 2giXB\Rt ■ £tI = ^SI^bBi cos 0++ , (12) 

Likewise, 

Wav = 2gflB\Rl ■ £>l\ = 2gjUBBi COS . (13) 



Thus, there can be a change in the precession frequency if 
0-+ + 0++. If we assume the magnitude of the magnetic field 
remains unchanged when it is reversed, i.e., - - Bq, 
all of the magnetic field components defined above may be 
rewritten in terms of these angles; 
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Because of the possibility of both changing the relative angle 
and magnitude of the magnetic field upon reversal of either 
the applied electric or magnetic fields, there are many poten- 
tial underlying causes for a non-zero derived value of a Bi^ak, 
Bcorr, or B„r field parameter Note, however, that the sim- 
plest way outside of an electronics issue to generate a non- 
zero Bcorr would be to have misaligned fields. 

B. N-odd combinations 

There are four other combinations that can be formed from 
the measured frequencies. These remaining combinations all 
contain the frequency difference between the A^ = +1 states, 
i.e., they are odd under reversal of A^. For a given state of the 
electric field and magnetic field, the difference of the preces- 
sion frequencies for the two A^ states will be given by equa- 
tions similar to those for CL)av.j,k contained in Table I, except 
2g will be replaced by Ag. In the limit where the molecule is 
fully polarized, which is very well satisfied under our condi- 
tions, Ag/g = k\S\ (see Eq. 2). 

Let us consider a simple model where the various states of 
the electric field are parametrized by three magnitudes. The 
first parameter will be the average electric field magnitude 
Sen, > 0. If there is a fixed charge density on the walls of the 
vapor cell, the electric field may change its magnitude from 
Sav when the direction is reversed. Let us denote this non- 
reversing component as £,„ . In addition to this non-reversing 
component of £, let us consider the possibility that the electiic 
field magnitude changes by an amount Scorr when the abso- 
lute magnetic field relative to the chosen laboratory frame is 
reversed". A full table of this parametrization of the mea- 
sured electric field in terms of these components for various 
states of the applied electric and magnetic field is shown in 
Table I. 

These three parameters that describe the electric field com- 
bine with the four parameters to describe the magnetic field 
to become the seven field parameters used in our experiment. 
Eight measured frequency combinations allow us to determine 



As with Bcorr, om cause of Scorr could be unwanted electrical coupling 
between the magnetic and electric field power supplies. 
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Label 



Physical Quantity 



A_ 
A, 

A+, 
A+, 
A_ 
A_ 



16gyUBS,„, 

SkBavgflsBay + U8„rglJ.BBcorr + ^k&corrg^i bB„, 

16gyUsB„,- 
l6gHB Bleat 

SkSavmsBnr + SkSnrgfJ-sBleok + ^kScorrgl^BBay 
SkSaygfiBBcon- + SkS„rgflBBay + SkScorrgf^BBleak 



A_ _ _ai ik&gyg^BBleak + MG,„.gflBB„,. + MScorrgfiBBcon " Idd^Seff 



TABLE n. Table of frequency combinations and associated field pa- 
rameters, along with the eEDM value dg. 



seven parameters describing experimental conditions, plus the 
value of the eEDM. This parameterization therefore repre- 
sents the maximum information that can be extracted from 
these frequency measurements alone. 

Using this parameterization for the four states of the mag- 
netic and electric fields, A_ + +w will be given by 



= k{&„t, + S„r + &corr)gf^B(Bav + Bcrr + Bleak + B,,,) 
+k{&av - £>nr + £>corr)glJ^B(Bav - Bcorr " Bleak + B,,, ) 
+k(Say + &nr " £>corr)gl^B(Bav + Bcorr - Bleak - Bnr) 
+k(Sav - ^nr " &corr)g^lB{Bav - Bcorr + Bleak — B„r) 

= 8kSavgl^BBav + SkS„rgfiBBcorr + SkScorrgl^BB„r -(15) 

In the above equation, each line corresponds to the applied 
electric and magnetic fields being in a different state (see Ta- 
ble I). In general, the misbehaving components of the mag- 
netic and electric fields are small compared to the average 
fields, i.e., \6„r\,\&corr\ « &av and \B„r\,\Bcorr\ « Bay Thus, 
Eq. 15 should be dominated by ^kSgvgf^BBav 

In a similar manner, A_,+,_w may be expressed as 

A-,+ -W = MSavg/^BBnr + SkSnrgflBBleak 

+U6corrgf^BBav , (16) 

and likewise A_,_,+w may be expressed as 

A- = MSavgf^BBcorr + SkSnrgHsBav 

+ SkScorrgf^BBleak ■ (17) 

This system of equations can be better expressed in terms of a 
matrix equation. 



ik 
2 



A+.+,+dLi A+._.+w A+,+ _w 
A+.+.-w A+_._(y A+,+,+(j; 
A+ _ +aj A+ + +w A+ _ -oj 















&corr 




A-,+ 




J 






.A_,_ 





(18) 

By solving this matrix equation, the electric field parameters 
can be determined from the measured frequency combina- 
tions. 

In addition to the electric field parameters, the eEDM also 
has an eff'ect on the measured precession frequencies. In par- 
ticular, the eEDM causes the frequency difference between the 

= +1 states to be shifted by -4sign(£)sign(B)(if,£e//. Thus, 



for the above frequency combinations, the eEDM signal can- 
cels. But for the frequency combination that is odd in A^, £, 
and B, we find 

A-,-,-W = Zk&avgUBBleak + SkSnrgfJ-BBur + 

UScorrgl^BBcorr ' IMe&eff ■ (19) 

A full listing of the frequency combinations is given in Ta- 
ble II. 

In addition to the eEDM signal, we find three additional 
terms that can produce a non-zero A_____w. The middle 
two terms, MS„rgfJ-BB„r and MScorrgl^aBcorr, generally cre- 
ate small contributions to A_ _._fa;, as they are products of a 
small component of the electric field and a small component 
of the magnetic field. Of these three additional terms, the term 
dikSavgl^sBieak IS generally anticipated to give the largest spu- 
rious contribution to A_ _ as it depends on the product of a 
misbehaving component of the magnetic field and the average 
magnitude of the electric field. However, compared to atomic 
experiments based on Eq. 1 (i.e., without the internal comag- 
netometer feature that arises from the n-doublet structure), 
the systematic error due to Bie„k is suppressed by kSav ~ 10"^. 

To demonstrate this systematic rejection, we manually cre- 
ated large misbehaving components of the magnetic and elec- 
tric fields and measured their impact on the eEDM signal 
channel. The results of this study are shown in Fig. 7. The 
prediction of Eq. 19 explains most of the correlation in the 
data (~ 98% for Bieak and ~ 80% for Bcorr and B„,); however, 
the fits in Fig. 7 show reduced values significantly larger 
than unity. The statistical uncertainties from the experiment 
are the only contribution used to compute these values, and 
therefore they might suggest suggest an additional systematic 
effect not accounted for in Eq. 19. 



V. THE EFFECT OF FIELD GRADIENTS 

While the above discussion shows that there is significant 
power in the internal comagnetometer to determine the pres- 
ence of systematic effects, the implicit assumption is that 
the applied fields are uniform and therefore identical for all 
molecules in the experiment. Here we extend our discussion 
to address the following question: to what extent do mag- 
netic and electric field gradients affect the field parameters of 
Sec. IV and A_ _ the eEDM signal channel? 

Perhaps the best proxy that exists for the size and strength 
of the field gradients in our experiment is the decay rate of the 
quantum beats F = l/T^, as defined in Eqs. 6-7. Two com- 
ponents contribute to the decay rate: the rate of dephasing 
of the beats due to all homogenous broadening effects I/T2 
(such as spontaneous emission and collisions) and the rate 
of dephasing of the beats due to field inhomogeneities l/r. 
The former rate is strictly speaking unknown; however, we 
take the fluorescence decay rate of the a state as a lower limit 
{I/T2 > 1/Tfl, where t„ x 50 fis; see Eq. 8)'^. For simplic- 
ity, we shall assume that F= I/T2 + 1/t. With larger field 



■ The decay rate of the fluorescence itself has several components, e.g. the 
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FIG. 7. (Color online) Shift in the eEDM channel, A_,_,_w as a func- 
tion of various products of field parameters, when such parameters 
are deliberately made large. Each plot corresponds to a different 
product of misbehaving components of the electric and magnetic 
field. The solid, blue lines correspond to a slope of one, which is 
the prediction of Eq. 19 based on the units of the x and y axes, with 
no adjustable parameters. 

gradients and therefore a larger distribution of precession fre- 
quencies, r becomes larger If we assume that the additional 
rate of dephasing beyond V - 1/ 7^2 is due to these gradients, 
the fractional field inhomogeneity can then be estimated from 
the width of the peak in the the Fourier spectrum. Namely, the 
width is given by 

6B 6lo 2 2 1 1 \ 

= — = = — r (20) 

Bay ojQ TWO m\ T2I 

where 58 is the full-width, half maximum of the distribution 
of magnetic fields and Bay is the average magnetic field, re- 
spectively, that is experienced by the ensemble of molecules. 
Likewise, 6a> = 2/t is the full-width, half maximum of the 



natural lifetime of the state and quenching of the electronic state due to 
collisions with the walls or other PbO, PbO*, Pb02, Pb202, etc. molecules. 
However, it is possible that some collisions can cause dephasing, without 
quenching the a state of the molecule. 
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FIG. 8. Example electric (left) and magnetic field (right) profiles 
given by Eqs. 21 and 22. For the electric field, a = 0.2 and for 
the magnetic field b = 0.2. Note that only the region < x < / is 
included because each detector in the experiment primarily captures 
fluorescence photons from one half of the cell. 

resulting distribution of precession frequencies and wo is the 
resulting average precession frequency. 

Equation 20 can be used to estimate an upper bound of the 
size of the gradients the ensemble of molecules experiences. 
Without an electric field applied, only the absolute inhomo- 
geneity of the B-field contributes to r, and a typical value of 
r = 1 /(37 jUs) for the beat decay rate at loq -2nx (0.75 MHz) 
is observed. Such a value, when combined with our upper 
Umit of T2, implies that 6B/Bav < 0.29%. At the much 
larger magnetic field = 2;7r x (2.25 MHz), the decay rate 
increases to F - 1 / (36 jUs), which actually implies a de- 
crease in the upper limit of 6B/Bav to < 0.11%'^. When a 
maximally homogeneous electric field that is large enough to 
polarize the molecules is applied, the typical decay rate in- 
creases from F ~ 1/(35 jUs) to F ~ 1/(30 yus) with an average 
Wo ~ 27r X (3 MHz). Because of Eq. 5, this implies that an 
application of an electric field increases the fractional inho- 
mogeneity in the value of B ■ £ from <0.09% to < 0.13%''*. 

Given that field gradients are observed in our experiment, 
the relevant question becomes to what extent do these gradi- 
ents impact the measurement of the physical quantities listed 
in Table II? To answer this question, let us first consider an il- 
luminating, yet simple model of a two-dimensional cross sec- 
tion of the gas of PbO, where electric field plates are posi- 
tioned parallel to the x axis to generate an electric field in the 
z direction. If the electric field plates are not infinite in ex- 
tent, the electric field will experience fringing. The effect of 



These results imply that T2 is smaller than our upper bound, or the B-field 
inhomogeneity is not due primarily to the applied field, or some combina- 
tion of the two. If T2 were equal to our upper bound and the field gradient 
was caused only by the applied magnetic field, then SBIB„r would be con- 
stant. 

Here "maximal homogeneity" corresponds to conditions where £ < and 
there is no electron emission as shown in Sec, II. This configuration is 
discussed in more detail below. 
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fringing will be to reduce the strength of the electric field as 
the edges of plates are approached. Assuming that the reduc- 
tion scales as to lowest order and the solution must obey 
Maxwell's equations, the resulting electric field profile will 
be given by 



£ = £(1 



1 



2boa—x . 



(21) 



where £o is the strength of the electric field at x - Q and 
z = 0, / is the size of the area in which this approximation 
holds, and a is a dimensionless parameter that describes the 
inhomogeneity of the field in the volume defined by /. Be- 
cause each detector in the experiment sees roughly half of the 
horizontal extent but the full vertical extent of the vapor cell, 
X and z must be constrained to < x < Z and -I < z < I- 
Let us also assume that there is an applied magnetic field in 
the z direction with a dB-Jdx gradient. To satisfy Maxwell's 
equations, the field must be given by 



x\ z 
■ bjjz + Bob-x , 



(22) 



where b describes the strength of the dB,/dx gradient and Bq 
is the magnetic field at x = and z = 0. These two field 
profiles are shown in Fig. 8. 

Neglecting the contributions from components of the mag- 
netic field perpendicular to the local direction of S as well as 
the eEDM terms, Eq. 5 yields simple expressions for the dif- 
ference and average frequencies of the = +1 states. The 
measured difference frequency of the = +1 states will be 
given by the average difference frequency over the volume, 
i.e.. Old = w+ - = (IgfisiB ■ S)kS) = Ikgfisi^ ■ £>) where 
the brackets denote averaging over all space. For the field gra- 
dients above, this averaging yields 

ojd = IkgfiB^ J '^^ J dzB- S 





^ u\ 1 




— b a 
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(23) 



Similarly, the average frequency of the = ± 1 states will be 
given by cjav = + cjJ) = Igjisi^ ■ &} or 



(24) 



The above integral can be approximated by expanding 1 /£ 
in a power series in the gradient a and then integrating. The 
result of the integration is 



._^_2(5H^ 
2 3 45 



(25) 



Thus we see that both the average frequency and the difference 
frequency can be affected by the presence of gradients. 

The determination of the various electric field parameters 
of Sec. IV is predicated upon the relationship between the dif- 
ference and average frequency, i.e., oj^ - k&ojav or a>+-a>- - 
kSi{LiJ+ + (jjJ) (see Eq. 2). Inserting Eqs. 23 and 25 into this 



relationship allows us to write an expression for the average 
derived electric field magnitude &d^,,- - to^/iktOav). For our 
example field configuration, the derived electric field is given 
by 



-*der 



1 
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6(2 + b) 
(40 -I- 4-4-b + 7b-)a^ 
45(2 + fc)2 



+ O(fl^) 



(26) 



Unhke what is assumed in Sec. IV, it is now apparent that the 
electric field magnitudes inferred from the data can depend 
on the strength of the background magnetic field gradient b. 
Moreover, the accuracy of the average electric field magnitude 
measurement also depends on the electric field gradient a. To 
demonstrate this, we calculate the actual average electric field 
magnitude. 
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(27) 



and take its difference with the derived electric field when 
b^O, 



->der ' 



704 
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Goa"^ + 0(a^) . 



(28) 



Therefore, with this field configuration, the measurement of 
the magnitude of the electric field is impacted by the electric 
field gradients at fourth order 

If the magnetic and/or electric field gradients change when 
the fields are reversed, various frequency combinations A,j 
can be non-zero, and therefore various field parameters such 
as Bcorr, Gcorr and Bleak will be inferred to be non-zero as well. 
In the context of our example, if we assume that the electric 
field gradients are unaffected by the magnetic field state and 
vice versa, we can define aj to be the electric field gradient 
parameter when the electric field has sign j. Likewise, let us 
define b^ to be the magnetic field gradient parameter when the 
magnetic field has sign k. It can then be shown that for the 
field profile used above 



bS) + a^{b+ - M] 
[(a+ - a_)b+ + (a+ - a_)i>_] 



leak 



Bq 
3 

--^ [a+b+ — a-b+ — a+b^ + a^b^] 



(29) 
(30) 
(31) 



Note that if the magnetic field gradients are reversed perfectly, 
i.e., b+ - b-, Gcorr = 0. Likewise, if the electric field gradi- 
ents are reversed perfectly, then = a+ and B^orr = 0. To 
make Bicak - 0, either the electric field gradients or the mag- 
netic field gradients must be reversed perfectly. Therefore, to 



" Note that for space considerations, tlie relevant term for this calculation 
is not shown in Eq. 26. 
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generate an apparent nonzero value of Bi^ak through field gra- 
dients, both field gradients must change upon reversal. 

One can imagine field profiles more complicated and more 
realistic than those used in the example above. Such mod- 
els could include misaligned fields (i.e. transverse field com- 
ponents), field gradients, and changing average magnitude of 
the fields. For example, using a misaligned magnetic field 
together with Maxwell's conformal solution for the fring- 
ing electric field in a parallel-plate capacitor [39], it can be 
shown that the derived value of the electric field magnitude 
will change as the magnetic field is rotated from the verti- 
cal direction. The exact quantitative relationship between the 
magnetic and electric field parameters of Sec. IV and the field 
gradients depends strongly on the exact nature of the chosen 
field profile in any given model. Because accurate field pro- 
files cannot be determined with this apparatus, construction of 
an accurate quantitative model for this experiment is impos- 
sible. Nevertheless, all of the models tested show the same 
qualitative behavior; changes in the magnetic field profile can 
aff'ect the derived electric field parameters and vice versa. To 
create a non-zero Scon- parameter, the magnetic field gradient 
must change when the magnetic field is reversed. Likewise, 
to create a non-zero Bcon- parameter, the electric field profile 
must change when it is reversed. Lastly, to create a non-zero 
Bleak parameter or a false eEDM signal, both the electric and 
magnetic field profiles must change when the respective fields 
are reversed. 

While the situation may now appear intractable, it remains 
important to note that if the gradients change upon the re- 
versal of the electric and/or magnetic fields, the beat decay 
rate F must also change. Therefore, forming combinations 
of F, analogous to the combinations of frequency discussed in 
Sec. IV, yields a quantitative measure of the amount of change 
in the field gradients. For example, the combination A+ +._F 
quantifies how much the magnetic field gradient (or, rather, 
the inhomogeneity in B ■ £) changes when the magnetic field 
is reversed. 

An example of this decay rate measurement and the sub- 
sequent combinations is shown in Fig. 9. As a function of 
an applied, fixed dB./dx gradient (henceforth, the term fixed 
shall refer to a gradient that does not reverse with its corre- 
sponding field, e.g., b+ = -b-), the average decay rate A+.+.+F 
changes and minimizes near zero gradient. However, the dif- 
ference in the decay rate when the magnetic field is reversed 
varies approximately linearly with an applied, fixed dB~/dx 
gradient, as shown in Fig. 9b. Note that this difference in the 
decay rate when the magnetic field is reversed does not go to 
zero when the average decay rate A+ +,+F is minimized, nor 
does it become zero for the two different detectors at the same 
applied dB^Jdx gradient'^. Lastly, a non-zero A+ _ _F is also 
observed and is shown to have a dependence on an applied, 
fixed dB./dx gradient. This indicates a change in decay rate 
upon reversal of both the magnetic and electric fields. 

In principle, one can apply not only a fixed gradient, but 
also a gradient that reverses with the magnetic field (hence- 
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' This implies the presence of another magnetic field gradient in the system. 
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FIG. 9. (Color online) Change in combinations of the beat decay rate 
as a function of an applied, fixed dB^/dx gradient. The red circles 
show the signal from the detector fixed in the +x direction and the 
blue squares show the signal from the detector along the -x direction, 
(a) The combination A+ +,+r is eight times the decay rate averaged 
over all state and field configurations. The combinations (b) A+ + _r 
and (c) A+,_,_r show how the decay rate changes with reversal of the 
magnetic field and both the magnetic and electric fields, respectively. 
The red, solid and blue, dashed lines are quadratic fits to the signals 
from the two detectors along +x and -x, respectively. 



forth, the term reversing shall refer to a gradient that does 
reverse with its corresponding field, e.g., b+ = b-). A 
deliberately-applied, reversing component may be useful in 
order to make the applied magnetic field more uniform. For 
example, if there is a displacement of the cell from the center 
of the Helmholtz coil, a linear field gradient that reverses with 
the applied field generated by the Helmholtz coil can shift the 
maximum of the field back to the center of the cell. Inside 
a magnetic shield, a Helmholtz coil will generate a field that 
has a quadratic gradient (e.g., a non-zero d~B-/dx^). If the 
geometric center of the cell and the center of the coil do not 
match, a linear gradient can be applied to eff'ectively shift the 
maximum of the field to the center of the cell. In doing so, 
such a linear gradient would have to be reversed with the field 
generated by the Helmholtz coils in order to keep the shifted 
maximum in the center of the cell. 

Because non-reversal of both electric and magnetic field 
gradients can shift the measured values of various parameters, 
a careful study was performed to determine the eff'ects of var- 
ious applied fixed and reversing magnetic field gradients on 
the field parameters of Sec. IV. By deliberately applying both 
fixed and reversing components of the magnetic field gradi- 
ents, we can quantitatively measure the effect of any partic- 
ular gradient on A+ _ _r, &corr^ Bi^ak, etc. Shown in the top 
panel of Fig. 10 is one example of the results if a magnetic 
field gradient is varied. Here, a dependence of &corr on a fixed 
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FIG. 10. (Color online) (a) Variation of &corr as a function of fixed 
dB-ldx. (b) The variation in eEDM signal channel A_ _ _toi as a func- 
tion of fixed dB-Jdz. The red circles show the signal from the de- 
tector fixed in the +x direction and the blue squares show the signal 
from the detector along the -x direction. 



dB-ldx gradient is observed. Such a dependence of &corr can 
be qualitatively explained using the toy model above, i.e., be- 
cause the gradient does not reverse with the magnetic field, we 
set - -b+ and Eq. 29 becomes - £o(fl+ + a-)b+/6. 

A particularly interesting quantity is A+ _ +r, which quan- 
tifies how the gradients of the electric field change upon its 
reversal. Thus, A+___+r can show the effects and general be- 
havior of electron emission described in Sec. II. Consider the 
data shown in Fig. 11. At the start of the run, we observe a 
large current due to electron emission that becomes smaller 
with time but trends to a non-zero, steady-state value of ap- 
proximately 1 .5 /iA. The time dependence of A+ _ +r follows 
a similar function form, but starts at a zero value and trends 
to a non-zero, steady-state value. This appearance of a non- 
zero A+ _ +r indicates that there is a distortion in the electric 
field profile that is correlated with electron emission reaching 
its steady-state'^. As the inhomogeneity of the electric field 
increases, a correlated change in Bcorr is observed. In the con- 
text of the gradient example described above, such a B^orr can 
arise if a+ + Such a B^-orr can also arise if electron emis- 
sion leads to changes in the angle of £, causing Q^^ + 9++ 
as described in Sec. TV. Some combination of these scenarios 
must be expected due to the non-reversing nature of the elec- 



While the cau.sal relationship is not cleai; it is possible that electron emis- 
sion can lead to an equilibrium state where various insulators in the cell 
have trapped charges on their surfaces or unwanted voltage drops due to 
cun'ent flow. Both of these effects could distort the electric field in such a 
way as to reduce the electric field neai' the emitting surface (thus reducing 
the emission current) but at the same time cause distortion of the electric 
field sunounding the molecules. 



-1 

-2 

< 

-4 
10 

B 

^ 5 

< 





o 



-50 



400 



^ h 



(b) 



H 1 1 1 1 1 

■ ; : : I (c) 



50 100 150 200 250 

Time since start of run (s) 

FIG. 11. (Color online) Change, as a function of time, in (a) the 
electron emission current, (b) A+ _ +r, and (c) S^o,, , when £ > 0. The 
electron emission cuixent is defined as the sum of all of the currents 
flowing into or out of the cell through the electrodes and is assumed 
to behave in a similar way to Fig. 4. At times f < 0, £ = 0. 



tron emission current, combined with the voltage drops and 
associated £-fields that accompany this current. 

Compared to magnetic field gradients, the electric field gra- 
dients cannot be as well controlled in this apparatus. However, 
it was empirically discovered that the size of the electric field 
gradients depends on the length of time the electric field is on 
during a data run (Te^o) compared to the length of time when 
the electric field is not applied {Te=o)- Reducing the duty cy- 
cle of the applied electric field, defined as T^^o /{Te=o + T^^o), 
below 30% reduces the size of A+ _ +r and the misbehaving 
field parameters £„,., Beon- and Bi^ak- 

The effects of both fixed and reversing components of 
dB,/dz, dB./dx, dB./dy, dB-ldx, dBy/dy, B^ and By on all 
the field parameters of Sec. IV and on the eEDM signal chan- 
nel A_,_,_w have been quantified. Two of the 14 components 
mentioned are observed to impact A_ _ The first is a fixed, 
transverse field B^, which our simple model above does not 
incorporate. The second and only magnetic field gradient is 
dB,/dz, as shown in Fig. 10b. While the simple model above 
also does not incorporate a dB./dz magnetic field gradient, the 
following linear gradients reproduce such an effect: 



B 



Boil+bj^z-bBoj^ 



z - a&ojx , 



(32) 
(33) 



where the range of x and z is specified by -/ < x < I and 
-I < z < I, respectively. With this configuration. 



cod - IgkfiB&oBo [1 + -ab 



(34) 
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Once again, if we assume the electric field profile is unaffected 
by the magnetic field and vice versa, then we can write a, (bj) 
as the size of the electric (magnetic) field gradient when the 
field has sign /. Then, 

4_ 

A_,_,_w - —gkfiB£>oBo(a+b+ - a-b+ — a+b- + a^b^) . (35) 

With a fixed dB^/dz, b- = -b+ and 
g 

- -gk^iB&QBQ{a+ - aJ)b+ . (36) 

Therefore, if the electric field gradient d&/dz does not re- 
verse perfectly, it is not surprising that a fixed dB-jdz gradient 
should produce a false eEDM signal. However, minimization 
of dB-jdz and other gradients can be achieved by minimizing 
the average beat decay rate A+ + +r, as described in the the 
next section. 



VI. EDM LIMIT 

Many types of experimental imperfections can conspire to- 
gether to create non-zero values of the parameters describing 
the magnetic and electric fields, as enumerated in Sec. IV. 
Parameters such as B^orn B„,. and £„^ can be generated by 
changing gradients, misalignment of the fields, non-reversing 
components of the magnetic or electric fields, or some com- 
bination thereof. However, parameters such as Scrr should 
only be generated by inhomogeneities in the magnetic or elec- 
tric field profile. Given that a systematic eEDM signal can be 
generated by gradients (e.g., by a fixed dB-jdz gradient, as 
shown in Sec. V), one must be careful to optimize the field 
profiles prior to taking any eEDM data. 

A. Minimizing gradients and optimizing ttie fields 

In general, the gradients of the magnetic field tend to be 
minimized when the overall lifetime is maximized. Using data 
similar to that shown in Fig. 9, we attempt to minimize the 
average beat decay rate A+,+,+r using both fixed and revers- 
ing components of applied magnetic field gradients. Given 
that we can apply both fixed and reversing components of five 
first-order gradients with our apparatus, this represents a diffi- 
cult optimization problem for which a solution cannot always 
be found. Moreover, the optimal solution may not make other 
combinations of the decay rate, such as A+ + _r or A+ _ +r, 
equal to zero. Such a situation indicates that although the gra- 
dients are minimized, they may not be identical upon reversal 
of the magnetic or electric fields. 

In order to gather useful eEDM data, the information 
gleaned above suggests a procedure for minimizing the spu- 
rious effects due to gradients. By applying both fixed and 
reversing components of both dB-Jdx and dB,/dy, we can 
demonstrate conditions that minimize the gradient (minimize 
A+ + +r) while equalizing the frequency measurements from 
the two detectors and driving A+,+,_r — > 0. Moreover, min- 
imizing A+,+,+r with both fixed and reversing components of 



dB-jdz ensures that dB-jdz is minimized in both the positive 
and negative field configurations. However, A+,+,+r is close 
to minimized with no applied dB-jdz gradient; therefore, the 
applied dB-jdz was set to zero in the course of the eEDM 
data set. These measures ensure the best magnetic field profile 
achievable. To control electron emission and hence minimize 
its contributions to non-reversing electric field components, 
the temperature of the vapor cell is maintained near 665 °C, 
the lowest temperature at which good signal to noise can be 
achieved. In addition, the electric field duty cycle is set be- 
low 30%, in order to ensure sufficient time for the non-Ohmic 
component of the leakage currents to settle to a small value. 
These two measures help to control the change in the electric 
field correlated with electron emission. 

Moreover, the derived values of the field parameters pro- 
vide a measure of how well the fields and their respective 
gradients reverse. In particular, a non-zero &corr or B„, can 
indicate a magnetic field profile that changes upon reversal. 
Likewise, a non-zero £„r and can indicate an electric 
field profile that changes upon reversal. A non-zero Bi^ak pa- 
rameter could indicate the non-reversal of both magnetic and 
electric field profiles, for if even one reversed perfectly, this 
term would be zero. 



B. Data collection 

For a given data run, defined as data taken under the same 
experimental conditions, the data is collected according to a 
relatively standard procedure. In order to calculate the field 
parameters of Sec. IV and the eEDM for a given data run, four 
applied field configurations must be used, corresponding to 
the four combinations of the signs of the magnetic and electric 
field. In general, an applied field configuration is selected, and 
512 laser shots (representing 5.12 s of data) are recorded with 
that field configuration. 

For every laser shot, each quantum beat signal is fit to ei- 
ther Eq. 6 or Eq. 7, depending on the data taking mode'^. The 
resulting collections of 512 best-fit values (e.g., F, (jl>\, ai and 
(^i) are then binned, with the number of bins determined by 
Scott's normal reference rule [40], and the resulting distribu- 
tion fit to a Gaussian. These distributions of best fit parameters 
show some outliers, which we believe are caused by fluctua- 
tions in the laser's output intensity'^. For each fit parame- 
ter, outliers are determined by Chauvenet's criterion [40] and 
these laser shots were excluded from the binning and averag- 
ing of all fit parameters. For a collection of 512 shots, approx- 



For data taken with microwaves, tlie N state selected for readout alternates 
from shot to shot, with the N state probed dining the first laser shot of 512 
chosen at random. 

For large fluctuations of the laser intensity, often the size of the beat signal 
will difi"er significantly from the average signal of the full 512 laser shots. 
The average signal is used to determine the initial guess of fit parameters 
(e.g., r, ai, a)\) in the non-linear fitting algorithm. For laser shots where 
the intensity drops significantly, the initial guess dift'ers from the best fit 
parameters, that it is not guaranteed that the non-linear fitting algorithm 
will converge to the best solution. 
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FIG. 12. (Color online) Final eEDM data set. (a) Best fit eEDM val- 
ues for each data run contained in the final data set. (b) The eEDM 
data shown as a function of electric field, (c) The eEDM data shown 
as a function of magnetic field. In each plot, the red circles represent 
data taken with two-beat technique, and the blue squares show data 
taken using microwave-erasure technique. The green, solid line in- 
dicates the best average A_.___tij, and the dashed, cyan lines indicate 
the 95% confidence interval for that average value. 



C. Statistics & data constraints 

Using data taken under the optimal conditions as a guide, 
we constrain our full collection of data runs in an at- 
tempt to eliminate any of the spurious effects that were de- 
scribed in Sec. IV-V. In order to ensure somewhat accu- 
rate reversal of the magnetic field gradients, we require that 
|A+,+,_r/A+,+,+r| < 0.0155 and \&corrl&m\ < 0.1%'°. To en- 
sure somewhat accurate reversal of the electric fields, we re- 
quire \&„rl&av\ < 0.1%. 

Using these constraints, a total of 4 hours of data is used 
in the final eEDM data set. This data includes two-beat data 
in the range of 0.65 G < Bav < 0.95 G (or, equivalently 
1.7 MHz < 2g^iBB„, < 2.4 MHz) and 100 V/cm < < 
125 V/cm. It also includes microwave data taken with similar 
fields. Lastly, some data where the procedure of Sec. VIA 
was either not applied (26% of the final data set) or paitially 
applied (20% of the final data set) manages to pass these cuts 
and therefore is included in the final data set. The final set of 
data that passes all cuts is shown in Fig. 12. 

The Ajj^itOi and A, j^tP for each run are then averaged to- 
gether, weighted by their respective errors. The result for the 
eEDM channel is A_,_,_w = 27r x (0.20 + 0.91) Hz. The xl 
value for the fit is 1.19 for 35 degrees of freedom; the proba- 
bility for a larger to occur is approximately 21%. 

As a final note, there are no detectable differences between 
the microwave data and two-beat data. However, the total 
amount of microwave data included is a factor of ~ 20 below 
that of the two-beat data, making the statistical error approxi- 
mately a factor ~ 4 times larger. 



D. Systematic errors 



imately 10 laser shots are typically excluded due to their best 
fit fli or fl2 values being significantly smaller than the mean. 

After 512 laser shots are recorded, the field configuration 
is changed. Within a given data set, each of the four field 
configurations is repeated approximately 16-32 times. The 
particular temporal order of the four applied field configura- 
tions has been found to not affect the final result. For example, 
with some runs, the electric field was reversed N times with 
the magnetic field positive or negative, the magnetic field was 
reversed and the magnetic shields degaussed, followed by an- 
other N reversals of the electric field. For other runs, the four 
magnetic and electric field configurations were cycled sequen- 
tially, without degaussing of the magnetic fields. No differ- 
ence between the two sets of runs is observed. 

With all field configurations measured N times, the com- 
binations A,j;j-w and A/j-^tF are then computed as a function 
of time using the four closest-spaced collections of 512 laser 
shots with the four required field configurations. Choosing the 
four collections most closely spaced in time minimizes the ef- 
fects of long term drifts in the magnetic and electric fields. 



Even with the cuts described above, the data shows clear 
evidence of non-ideal electric and magnetic field reversals. 
Therefore, an estimate of the error caused by the presence of 
imperfections of the magnetic and electric field profiles must 
be made. There are two magnetic field imperfections, dB-Jdz 
and Z?i, that when coupled with electric field imperfections 
impact A_,_,_w substantially. In the case of dB^/dz, Eq. 36 
acts as a guide and suggests that, in the presence of a fixed 
magnetic field imperfection and an electric field gradient that 
changes with reversal of the electric field, the data should obey 

where SS is a measure of the changing electric field gradient 
and c', is the constant correlation coefficient. To mea- 
sure c'. we deliberately apply large values of dB,/dz and 
increase the electric field duty cycle to amplify the negative 



These particular values were chosen to include the most data while min- 
imizing the scatter in the field parameters and A_ _ For example, a 
change from |A+ + -F/A-,- + -i-FI < 0.015 to 0.016 increases the^J of a Gaus- 
sian fit to the distribution of A_____w by almost a factor of 2. 
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FIG. 13. (Color online) Observed correlation between A_____6o» and 
applied (a) and (b) dB^/dz while varying the electric field duty 
cycle, which impacts A+._.+r and therefore the reversibility of the 
electric field gradients. The green lines represent the best linear 
least-squares fit of Eqs. 38-39, and are used to extract the c--^edm 
and Cy EOM coefficients. In both(a) and (b), the red circles show the 
signal from the detector fixed in the +x direction and the blue squares 
show the signal from the detector along the -x direction. 



effect of electron emission on the quality of the electric field. 
Because there is no direct measurement of the imperfection of 
the electric field gradient reversal 5&, we use as a proxy the 
quantity A+ _ +r, which quantifies how much the beat decay 
rate changes when the electric field is reversed. We use a sim- 
ple first order approximation that the two are linearly propor- 
tional, i.e. A+__ +r « ^5&, where ^ is the first-order expansion 
constant. Figure 13b shows this data and the subsequent fit to 



dz 



(38) 



where c^^~^edm = C. 



ZZ.EDM 



l( is the tunable constant extracted 
from the fit and dB^Jdz, A+ _ +r and A_ _ _w are the mea- 
sured quantities. The fit yields a value of c~-^edm / C^tt) - 
160 ± 15 Hz-yL(s/(/iG/cm). We note in passing that the sim- 
pler relation A_ _ _w = Czz,edm^, which ignores the impact 
of the changing electric field gradient, provides a poor fit to to 
data as ^ is varied. 

A similar dependence is observed in the case of B^. Empir- 
ically, the dependence of A_ on is given by 



(39) 



where c^ edm is the correlation coefficient for a fixed B^ com- 
ponent to the magnetic field. This dependence is shown in 
Fig. 13a. The constant c^ edm is measured in an analogous 
way to c~.^edm, and yields Cx,edmI{'^^) - 391+71 Hz-yL(s/mG. 



To estimate the systematic effect due to the term on 

the eEDM data set, we first calculate a limit on the maximum 
size of any remnant dB-jdz during the eEDM measurement. 
Two methods were used to measure this remnant dB^Jdz- 
First, as described in Sec. VIA, we attempted to cancel any 
fixed, remnant dBJdz by applying a dB^/dz gradient to can- 
cel the remnant gradient and therefore minimize A+ + +r. This 
procedure found that the background dB./dz was close to zero 
and therefore we applied no dB-jdz gradient. Over the course 
of the eEDM data set, this measurement was repeated at var- 
ious stages and all measurements are included in computing 
the final, average value of dB-Jdz- The second method in- 
volves masking the excitation laser such that only molecules 
on the top half of the vapor cell are excited into the a state. On 
the subsequent 512 laser shots, the mask is switched such that 
only molecules on the bottom half of the vapor cell are excited 
into the a state. The applied dB-jdz that makes the difference 
in the measured frequency for molecules on the top and bot- 
tom halves of the cell zero is the applied dB-jdz that cancels 
the background magnetic field gradient. All these measure- 
ments are combined and yield an average fixed, remnant value 
of dB-Jdz = -7.6 + 8.1 yuG/cm. 

An estimate of the background, fixed component of B^ must 
also be obtained in order to calculate the systematic effect of 
this transverse field. Such a limit is taken from measurements 
of the misalignment of the electric and magnetic fields (see 
Sec. III). These measurements suggest that the background, 
fixed Bx field component is B^ - -0.38 ± 0.16 mG. Because 
this value for the background, fixed B ^ is consistent with zero 
at the 3cr level, no compensation was applied to eliminate it. 

For each run indexed by k, the measured value of A+ _ +r 
and the estimates of dB-jdz and B^ are used to determine their 
effect on the eEDM channel for that run, (5A- - -w)i (using 
Eqs. 38-39). To determine the mean effect of the systematic 
on the full data set, the {6l^---<jS)k are averaged together ac- 
cording to 



6 A w = 



2 



A ii),k 



z 



IS.---<i},k J 



(40) 



where cta ^.i- is the statistical uncertainty in the eEDM value 

for run k. Note that the final systematic error is determined by 
weighting the magnitude of the individual systematic errors 
by the statistical error for that run. Each value of {6l^---(jS)k 

has a corresponding uncertainty ct^a and using simple 

error propagation, we compute the resulting uncertainty in 
6h w as 



z 

V k 



SA co.k 



Z 



(41) 



Alternatively, the standard deviation of the weighted average 
can also be computed. Such a standard deviation can be used 
to estimate the uncertainty in the mean value. If there are N 
runs, this estimate of the uncertainty of the mean systematic 
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Source 


Formula 


Values 


Shift 


Error 


Fixed dB-Jdz 
Fixed S V 

Product of £„,, and 
Product of £„, and B„, 

Product of &corr and Bcon 


Czz.EDM i 
^kGcorrg^lsBcorr 


^ = -7.6 ± 8.1 yuG/cm 
= -0.38 ± 0.16 mG 
Gay varies 
£„,- = (51 ±9)x 10-3 V/cm 
= (5 ± 9) X 10-3 V/cm 


A+._.+r = (0.23 ± 0.20) 1/ms 
A+___+r = (0.23 ± 0.20) 1/ms 
Sm. = -0.01 ± 0.52 fiG 
B„r = 5.7 ± 8.0 nG 
Bcoir = -2.5 ± 2.4 yuG 


-0.18 Hz 
-0.034 Hz 
-0.004 Hz 
0.00002 Hz 
-0.00014 Hz 


0.16 Hz 
0.030 Hz 
0.048 Hz 
0.00046 Hz 
0.00011 Hz 


Total 








-0.21 Hz 


0.16 Hz 



TABLE 111. Systematic eiTor budget. Shown are contributions from each source considered in the analysis to A_,_,_6oi/(27r), in units of Hz. The 
sources shown are described in the text. 



effect would be given by 



6^ a 



N 



K k 



2 

I A- 



A a),k ) 

(42) 



For any given systematic effect, we take the larger of cr- 



■ as the final systematic uncertainty. 



A breakdown of the most important individual contribu- 
tions to the systematic error is shown in Table III. For errors 
caused by magnetic field imperfections, we explicitly include 
contributions only from the two largest observed sources of 
correlation with the eEDM channel, namely those due to a 
fixed dBJdz gradient or a fixed field. Once again, our anal- 
ysis indicates that these enter the eEDM channel due to their 
coupling to uncontrolled, non-reversing electric field gradi- 
ents. For reference. Table III also contains the systematic con- 
tributions to A_ _ _(y described in Sec. IV. These contributions 
are calculated directly from the data, namely by constructing 
all possible combinations Aiji-o). The corresponding physi- 
cal field parameters for each run, including Bieak, B,,,-, Scon-, 
etc., are computed from these combinations via Table II and 
Eq. 18. With the field parameters, the systematic error for 
each run is computed using Eq. 19. Eqs. 40-42 are used to 
average together all the field parameters and systematic errors 
for the full data set, and the results are shown in Table III. 

It is important to note that the largest systematic error from 
the type of uniform field parameters discussed in Sec. IV is 
%g^Bk&a\ Bieak- While Bii,ak Can be created by changing mag- 
netic and electric field gradients, it can also be generated sim- 
ply by a leakage current. Given that the measured leakage 
currents in the cell are of the order of 10 yuA, it is certainly 
plausible that these leakage currents could generate a B/paA 
within our error of 0.53 jjG. However, if this experiment 
was performed with a system with no internal comagnetome- 
ter, i.e., where the Larmor precession frequency is given by 
Eq. 1, the contribution to the uncertainty of the systematic 
shift from such a large leakage current would be a factor of 
1 l{k&a\) ~ 100 larger than it is here. Hence, without the inter- 
nal comagnetometer, the leakage current contribution would 
be the largest single contributor to the systematic error 

Finally, we note in passing that we also considered other 
possible sources of systematic error One example is possi- 
ble errors arising from differences between the density and/or 
velocity of the populations in the = +1 states. By chang- 
ing the detuning of the laser from the center of the Doppler- 



broadened line, we could selectively populate one state 
more than another By removing the retro-reflecting mirror 
and setting the laser detuning in the middle of the Doppler 
broadened line, each state would correspond to an equal 
and opposite velocity class, which can lead to a spatial sep- 
aration of the two populations. We found no dependence of 
A_ _ _w on these effects. 

We therefore quote the final values of A_ _ _w = 0.20 + 
0.22,y„+0.91stat+0.17sys,)Hzandt/, = (-4.4+9;5„at+1.8syst)x 
10"^^ e-cm, where stat denotes the statistical Icr error and syst 
denotes the systematic shift and its corresponding 1 cr error A 
limit on the magnitude of the eEDM is obtained by integrating 
the assumed underlying Gaussian distribution symmetrically 
about the mean value, with the standard deviation taken as the 
quadrature sum of the statistical and systematic errors. The 
result is \de\ < 1.7 x 10"^^ e-cm, at 90% confidence. 



VII. CONCLUSION 

We have demonstrated that the fl^E^(v = 5, / - \) Q.- 
doublet state of PbO is capable of good systematic rejection 
and control in an electron electric dipole moment search. Us- 
ing PbO, we have obtained a limit of \de\ < 1.7 x 10""* e-cm 
(90% confidence), only about a factor of 20 worse than the 
world's best experimental limit [10]. Given the presence of 
significant non-reversing electric and magnetic field gradients 
(the former apparently due to uncontrolled electron emission 
from the top electrode) and leakage currents on the order of 
10 juA, obtaining such a limit provides a clear example of the 
power and flexibility an O-doublet state gives toward diag- 
nosing and controlling systematic errors in this type of exper- 
iment. 

The primary reason for the demonstrated level of system- 
atic rejection stems from the use of the two states that allow 
for measurement of the average magnetic field in all config- 
urations of the experiment. For this reason, this Q-doublet 
structure has been referred to as an "internal comagnetome- 
ter". Compared to a traditional comagnetometer [41], this sys- 
tem is sensitive to exactly the same magnetic field with which 
the molecules used to detect the eEDM are interacting. More- 
over, the ability to accurately measure the electric field using 
the molecules, because of the dependence of the g factor on 
electric field, gives even more information than a traditional 
comagnetometer 

Perhaps the only difficulty with this type of level structure 
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Evolution of 
spin 




FIG. 14. Detection geometry for an eEDM experiment that measures 
phase. See Appendix A for description. 



is the breakdown of the comagnetometer function in the pres- 
ence of magnetic and electric field gradients, rather than sim- 
ply uniform fields. This breakdown is caused by the ensemble 
of molecules being distributed in a finite volume with non- 
uniform magnetic and electric fields. This complication can 
have an impact on the inferred field parameters such as the 
average electric field and the component of the magnetic field 
that mimics a field due to leakage currents. However, by uti- 
lizing all available information, such as changes in the quan- 
tum beat decay rate and the derived electric and magnetic field 
parameters, large effects on the eEDM signal can be avoided. 

The next generation of experiments that use molecules with 
similar level structure, such as those based on ThO [17], 
HfF+[18, 19] and WC [16], should, by extension, have sim- 
ilarly good systematic rejection. In these experiments, the 
breakdown of the magnetometer due to gradients should 
pose less of a problem, as none require the use of a high- 
temperature vapor cell and thus avoid the complications inher- 
ent to the experimental apparatus described here. For exam- 
ple, in the ongoing ThO experiment [42], which uses a molec- 
ular beam, the electric and magnetic fields are much more uni- 
form than in this work. Moreover, the fields are more well un- 
derstood in that experiment; therefore more realistic and de- 
tailed modeling of the shifts encountered due to gradients can 
be undertaken. Given the result presented in this paper, we ex- 
pect that these future experiments should obtain dramatically 
better systematic rejection than that obtained here. 
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Appendix A: Systematic rejection when measuring phase 

In case of an experiment detecting phase precession in a 
beam (e.g., the ThO eEDM experiment [17]), most of the 
discussion of Sec. IV still applies; however, special atten- 
tion must be dedicated to how the phase angle is detected. 
In some experiments, it is possible to determine not only the 
total amount of phase precession, but also the direction of the 
rotation. This is the case in Ramsey's method of separated 
oscillatory fields (SOF) when one uses rotating fields at the 
start and end of the sequence. For this discussion, however, 
consider the case where the direction of spins' precession is 
detected by rotating the detection coordinate system relative 
to the preparation coordinate system by an angle 9. This is 
possible in the context of Ramsey's method of SOF if the two 
oscillating fields have different phases. 

In particular, consider an experiment where one detects the 
phase angle in a coordinate system defined by the axes x' and 
y that is rotated by an angle 6 with respect to the coordinate 
frame where the x coordinate is defined by the initial orienta- 
tion of the spin, as shown in Fig. 14. After evolving for a time 
T, the spin will have precessed by an angle (p+ if the spin is 
rotating in the clockwise direction. Let us assume that clock- 
wise rotation corresponds to Z? > 0. If the B field is negative, 
the spin will therefore rotate through an angle 4>- - -<P+- Af- 
ter precessing for a time t in a negative (positive) B field, the 
measured angle in the detection coordinate system will be a_ 
{a+). In this case, it is straightforward to show from Fig. 14 
and the parameterization used in Table I that the various aav,j,k 
are given by 



= 2gHB{B„, + Bcorr + B/eak + B„,)t - (Al) 

= 2g^iB{B„, - B,,,.,- - Bu,„k + B„,)t - 9 (A2) 

= 2glUB(Bav + B,orr - Bleak " B„,)t + d (A3) 

= IgUBiBav - Bcorr + B,eak - B„,.)t + 9 . (A4) 



The A^-even combinations then become 

A+_+_+a = l6gflBBavT 
A+-,+a = l6gflBBconT 

A+,+ -a = l6gfiBB„rT - I 

A+ - -a = l6gl^BBleakT . 



(A5) 
(A6) 
(A7) 
(A8) 



Note that adding an offset in the detection coordinates looks as 
if one is inducing a non-reversing component of the magnetic 
field. 

One possible way to detect the angle a involves projecting 
the spin along either the x' axis (with probability cos^ a) and 
the y' axis (with probability sin^ or) [17] . One then defines the 
asymmetry 



A = 



N, - Ny 
+ Ny 



cos^ a - sin^ a 
cos^ a + sin^ a 



= cos(2a) . 



(A9) 



Strictly speaking the A's are not additive and thus cannot be 
used to compute the combinations except in an approximation 
where a x jt/4-. In that case. 



(AlO) 
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Type Label Physical Quantity 



Phase A+_+ 


+a 


lei/JijBnvT 


A-.+ 


+« 




A+.- 


+« 




A+.+ 


_Q' 


l6-gtiBB„,T - W 


A+,- 




l6gflBBleakT 


A-,+ 




Sk6a„gflBB„rT + ^k8„,.gHBBleakT + ikScorrgl^BBavT 


A_._ 


+ « 


USaygl^BBcorrT + U6„,-glJBBavT + USn„,gflBBleakT 


A_._ 




%kS„„gtlBBl,„kT + SkG„rgllBB„rT + SkGcorrgf^BBcorrT - l6d,GeffT 


Asymmetry A+ + 


+ 


-32giiBSBT 


A-.+ 


+A 


-IGkBayg^BBayT - l6kS„rgflBBcorrT - \6k&,;„.rgHBB i„.T 


A+.- 


+A 


-32gHBBconT 


A+.+ 


_A 


-32gHBB,„T + 166 


A+.- 


_A 


-32gfiBBleakT 


A-.+ 


_A 


-\()k&„ygHBB„,T - l6kS„i-gllBBleakT- IGkB^onmsB cn-T 


A_._ 


+A 


-IGkBavmBBcorrT - l6k8„rgJJBBavT ~ IGkScongt^BBleakT 


A_._ 


_A 


-l6k8„yg^iBBi,„tT- l6k8i,rg^iBB„,T- \6kE,,„,,g^BBco,,T + 32d,&,ffT 



TABLE IV. Full phase and asymmetry combinations for an eEDM experiment that measures total precession through a magnetic field. 



and the average asymmetries become 



= -AgflBiBav + Bcorr + Bleak + B„,)t + 20 + - (All) 

Aav,-,+ = -4glJB(B„, - Baorr ' B^ak + B„r)T + 29+^ (A12) 

Aav,^,- = -^t^BiBar + B,„,, - B,eak - B„)t -20+^ (A13) 

Aav,-.- = -4glJB{B„, - B,orr + B,eak - B„,)t -20+^ (A14) 

Using the above expressions to form the A^-even combinations 
yields the following result: 

A+,+,+A = -32gHBBa.T + 4jt (A15) 

A+,_,+A = -32gHBBcorrT (A16) 

A+.+._A = -32gfiBB„rT + 160 (A17) 

K-,-A = -32gfiBB,eakT . (A18) 



One can then define 6B = Bay — {7i/4)/{2giJBT), which is the 
magnetic field magnitude that moves the average spin preces- 
sion away from a = n/4. With this definition. 



A+,+,+A = -32g^iB6BT , 



(A19) 



which more closely resembles the expressions for A+ + +(y and 

A+.+.+a. 

For the A^-odd combinations, the resulting expressions are 
almost identical to the case of measuring frequency. Consider 
first the case of detecting the asymmetry A; both the and 7t/2 
terms cancel in the A^-odd combinations, leaving expressions 
identical to those for A- y^w except they are multiplied by -2t 
Likewise, in the case of detecting the phase a, the constant 6 
terms cancel in the A^-odd combinations. Thus, the expres- 
sions for A_j,aQ' are identical to those for A_ except they 
are multiplied by t. All of these results are enumerated for 
reference in Table IV. 
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